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We study dark matter (DM) which is cosmologically long-lived because of standard model (SM) 
symmetries. In these models an approximate stabilizing symmetry emerges accidentally, in analogy 
with baryon and lepton number in the renormalizable SM. Adopting an effective theory approach, 
we classify DM models according to representations of SU{3)c x SU(2)l x U{1)y x 17(1)s x U{1)l, 
allowing for all operators permitted by symmetry, with weak scale DM and a cutoff at or below 
the Planck scale. We identify representations containing a neutral long-lived state, thus excluding 
dimension four and five operators that mediate dangerously prompt DM decay into SM particles. 
The DM relic abundance is obtained via thermal freeze-out or, since effectively stable DM often car¬ 
ries baryon or lepton number, asymmetry sharing through the very operators that induce eventual 
DM decay. We also incorporate baryon and lepton number violation with a spurion that parame¬ 
terizes hard breaking by arbitrary units. However, since proton stability precludes certain spurions, 
a residual symmetry persists, maintaining the cosmological stability of certain DM representations. 
Finally, we survey the phenomenology of effectively stable DM as manifested in probes of direct 
detection, indirect detection, and proton decay. 

I. INTRODUCTION 

Dark matter (DM) is elegantly accounted for by a neu¬ 
tral, cosmologically long-lived particle beyond the stan¬ 
dard model (SM). In most circumstances, however, DM 
stability is ensured by a symmetry that is simply im¬ 
posed by fiat. While as much can be expected from 
DM effective theories, similarly ad hoc choices are of¬ 
ten needed for their ultraviolet completions, for instance 
in theories of supersymmetry or extra dimensions where 
i?-parity or KK-paiity are assumed. 

In contrast, the SM implements stability with less con¬ 
trivance: charge stabilizes the electron, while angular 
momentum stabilizes the neutrino. The proton is not 
guaranteed to be stable, but like DM it is cosmologically 
long-lived, with current limits bounding its lifetime to 
be greater than ^ 10^^ years. Famously, the SM gauge 
symmetry explicitly forbids baryon and lepton number 
violation at the renormalizable level, suggesting an ar¬ 
gument for proton stability from effective theory. This 
mechanism, whereby an approximate symmetry arises 
as the byproduct of existing symmetries, is sometimes 
referred to as an accidental symmetry. 

In this paper we argue that DM, like the proton, can be 
cosmologically stable as an accident of SM symmetries. 

For our analysis we consider the symmetry group 

SU{S)c X SUi2)L X U{1)y x U(1)b x C/(1)l, (1) 

which is exactly preserved in the renormalizable SM to 
all orders in perturbation theory. Beyond the renormal¬ 
izable level, B and L may be approximate or exact, de¬ 
pending on whether they are gauged in the ultraviolet. 
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While SM quark and lepton flavor are approximate sym¬ 
metries, we will not consider them here. 

In our analysis, we enumerate models according to the 
quantum numbers of DM under the SM symmetry group. 
We take the stance of effective theory throughout: all 
interactions, renormalizable and non-renormalizable, are 
to be included with order one coefficients. We assume a 
DM mass M near the weak scale and an effective theory 
cutoff A at or below the Planck scale. 

To start, we discard all representations without a com¬ 
ponent neutral under color and electromagnetism. We 
then discard all representations that permit DM decay 
into SM particles via dimension four or five operators. 
In such cases even a Planck scale cutoff is insufficient to 
prevent cosmologically prompt DM decay. On the other 
hand, decays induced by dimension six operators hover 
right at the boundary of current bounds from indirect 
detection, assuming a cutoff near the scale of grand uni¬ 
fication (GUT) [T]. 

Applying these criteria, we enumerate models of effec¬ 
tively stable DM, focusing on all possible fermionic and 
scalar DM candidates whose leading decays occur at di¬ 
mension six or seven. These results are summarized in 
Tabs. HI and HIl for dimension six and Tabs. Hill and IIVI for 
dimension seven decays, respectively. 

Our approach resonates with that of minimal DM HE] 
but with the crucial difference that we incorporate both 
B and L and faithfully assume the presence of all inter¬ 
actions not forbidden by symmetry. Without B and L, 
representations smaller than the quintuplet of SU(2)l, 
DM will promptly decay to SM particles via renormaliz¬ 
able interactions US]. This is consistent with our own 
finding that effectively stable DM in small representa¬ 
tions of SU{2)l must carry B or L. For exactly this 
reason many of these models may be generated through 
the mechanism of asymmetric DM ESI. Previous au¬ 
thors have also considered accidental stabilization of DM 
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by flavor symmetries [7] or new gauge symmetries [5]. 

Finally, we consider the possibility that B and L 
are merely approximate. For this analysis we intro¬ 
duce a spurion parameterizing the hard breaking of 
baryon and lepton number by arbitrary units Ai? and 
AL, respectively. In many models of effectively sta¬ 
ble DM, this induces prompt decays. However, not all 
values of {AB,AL) are phenomenologically safe: the 
non-observation of proton decay suggests that dimen¬ 
sion six operators of the form Q^L should be forbidden, 
thus precluding hard breaking by units of {AB, AL) = 
(1,1). Remarkably, given arbitrary breaking by any 
{AB,AL) ^ (1)1)) we still maintain a handful of vi¬ 
able candidates for effectively stable DM. In these mod¬ 
els, DM is long-lived because of the SM gauge symmetry 
together with the stability of the proton. 

The remainder of this paper is as follows. In Sec. |lT] 
we enumerate viable models of effectively stable DM us¬ 
ing the criteria described above. In Sec. [ITT] we study 
experimental constraints on these theories from direct 
detection, indirect detection, and proton decay. Finally, 
we summarize our results and discuss future directions 
in Sec. |IVl 


II. CLASSIFICATION OF MODELS 


In this section we enumerate representations of the 
SM symmetry group with a DM component that is neu¬ 
tral, cosmologically long-lived, and generated with the 
observed relic abundance. We adopt an effective theory 
perspective in which any operator allowed by symmetry 
is present with a strength set by a sub-Planckian cutoff. 


A. Neutrality 


The quantum numbers of DM are parameterized by 
a discrete choice of representations for 51/(3)c and 
SU{2)l together with a continuous choice of charges for 
U{1)y, U{1)b, and For simplicity, we focus on 

pure gauge eigenstates, assuming that DM is either a 
complex scalar or a Dirac fermion. 

To begin, we restrict to representations that have a 
neutral component under the unbroken SM gauge sym¬ 
metry. Thus, DM is an SU{3)c singlet. For DM that 
is a fc-plet under SU{2)l, we require that k > 2|y| -|- 1, 
where the hypercharge Y is quantized to an integer or 
half-integer value if k is odd or even, respectively. 

Charge neutrality places no constraints on the B or L 
charges of DM. While irrational values of B and L are al¬ 
lowed a priori, this is literally equivalent to enforcing DM 
number as an exact symmetry of the Lagrangian. Fur¬ 
thermore, exact global symmetries are known to conflict 
with black hole no-hair theorems [H|. For these reasons 
we restrict to rational values B and L. 


While B and L are symmetries of the renormalizable 
SM, they may of course be spontaneously or explicitly 
broken in the full theory. To parameterize these ef¬ 
fects conservatively, we introduce effective hard break¬ 
ing of B and L into the low energy theory with a di¬ 
mensionless spurion for symmetry breaking by units of 
(AB, AL). For remainder of this section we assume that 
(AB,AL) = (0,0), but return to the issue of explicit 
breaking later on in Sec. |III| 


B. Stability 


We now determine the leading operators that are al¬ 
lowed by symmetry and mediate DM decay. An operator 
that induces DM decay takes the form 

Odm = AOsm ) (2) 

where X is the fermion or scalar held that contains DM 
and C>sM is an operator composed entirely of SM fields. 
For later convenience, we define 


N = [Odm] ) 


(3) 


to be the dimension of the DM decay operator. The 
quantum numbers of X are equal and opposite to those 
of OsM) so to enumerate all decay operators it suffices 
to determine all operators Osm of a given charge and 
operator dimension. 

We define a fiducial decay rate into SM particles. 


T{X ^ SM) 


47r \ A / 


(4) 


corresponding to two-body decay via Odm- Depending 
on the precise operator, this decay may be three-body or 
higher. Moreover, decays into SM fermions will involve 
flavor structures that may further suppress the width. 
In any case, the fiducial decay rate in Eq. 0 should be 
taken as an overestimate. 

By definition, a cosmologically stable DM particle X 
has a lifetime of order the age of the universe, 

t{X —)■ SM) > 10^® sec (age of universe). (5) 


However, this bound is weaker than experimental con¬ 
straints on cosmic ray production from DM decay into 
positrons lin!) gamma rays m, antiprotons m. and 
neutrinos m- These limits all place similar constraints 
on the DM lifetime, of order 

t{X —> SM) > 10^® sec (indirect detection), (6) 

for M of order the weak scale. While CMB bounds are 
also stronger the one from the age of the universe, they 
are still weaker than indirect search bounds by ~ 2 — 3 
orders of magnitude m- 
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Comparing Eq. ^ to Eq. Q , we obtain upper bounds 
on the cutoff of higher dimension operators. For di¬ 
mension five, six, and seven operators, this implies a 
schematic lower bound on the cutoff. 


A> 


M 

1 TeV 
M 

HhV 

M 

IlhV 


10^9 GeV , N = 5 
IQi® GeV, N = 6 
10^2 GeV, N = 7 


(7) 


Since Eq. is an overestimate, this bound on the cutoff 
is conservative. 

For dimension five operators, a sub-Planckian cutoff 
A would require M < 10 keV for DM, so hereafter we 
consider only DM decay via dimension six and seven 
operators. Notably, dimension six operators are a par¬ 
ticularly intriguing portal for weak scale DM because 
GUT-suppressed dimension six operators induce decays 
that lie just at the boundary of current experimental 
limits. For dimension seven operators, the bound on A 
is even smaller. Fig. shows the lower bound on A as a 
function of M for dimension six and seven decay oper¬ 
ators. The bounds for dimension five operators are not 
shown because they require a cutoff above the Planck 
scale. Also shown are the natural scale for DM mass M 
near the weak scale and cutoff scale A at the GUT scale, 
as well as regions excluded by experimental limits from 
LEP and direct detection, to be discussed later. 

To summarize, a sub-Planckian cutoff implies that cos- 
mologically stable, weak scale DM forbids all dimension 
four or five decay operators. This is a stringent con¬ 
straint on the DM representation. For example, since 
a pure singlet scalar X can couple to any dimension 
four operator C>sm in the SM, the associated DM will 
promptly decay at dimension five. At the opposite ex¬ 
treme is scalar or fermionic DM with extremely large 
charges. In this case lower dimensional operators are 
forbidden simply because so many SM fields have to be 
included in the decay operator just to preserve the sym¬ 
metry. As noted in [2], this occurs for fc-plets of SU{2)l 
with large values of k. 

All of our fermionic DM candidates and most of our 
scalar DM candidates carry non-zero B or L. The reason 
is that most representations with zero B or L typically 
have very low dimension decay operators. For example 
whenever Osm is a SM fermion bilinear that is B and L 
neutral, it is always possible to construct a lower dimen¬ 
sion operator via replacements of the form QU —>■ H\ 
LE —>• H, Qa^Q —> W^, etc. As a result, viable DM 

tends to carry B or L, or reside in a large representation 
of SU{2)l so that gauge invariance requires many Higgs 
fields in the operator. 

Large hypercharge can serve the same role in terms of 
stability, but as noted earlier, a neutral DM component 
requires that k > 2\Y\ -|- 1 for a fc-plet under SU{2)l. 
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Figure 1. Lower bound on the effective theory cutoff A as 
a function of the DM mass M in order to avoid cosmologi- 
cally prompt decay via dimension six (gold dashed) and seven 
(teal dot-dashed) operators. LEP (red shaded) excludes light, 
weakly interacting or hypercharged DM, while direct detec¬ 
tion (blue shaded) excludes hypercharged DM that cannot 
accommodate an inelastic splitting to evade bounds. 
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This in turn bounds the net SU{2)l x U{1)y charges of 
the SM fields that couple to X. In particular, it largely 
excludes decay operators involving E, since this SM field 
has a sizable hypercharge and no SU (2 )l charge. 

Tabs. [l] and [iT] list all fermionic and scalar DM rep¬ 
resentations, respectively, whose leading decay is medi¬ 
ated by a dimension six operator. Every representation 
is color neutral and every representation carries B or L 
except the scalar SU{2)l sextet. For the fermionic DM 
we forbid both members of the Dirac pair from decay¬ 
ing via dimension five or lower operators. Also, shown 
are various attributes of the model regarding direct de¬ 
tection and explicit B and L breaking, to be discussed 
later. Tabs. m and CYl are the same as Tabs. Uand El 
except they apply to DM representations whose leading 
decay is mediated by a dimension seven operator. 


C. Relic Abundance 

Let us comment briefly on the origin of the DM relic 
abundance. DM with non-trivial SM gauge charges will 
be equilibrated with the thermal plasma in the early uni¬ 
verse, in which case DM annihilations are mediated by 
gauge interactions whose strength is fixed by the charges 
of X. In general there may be additional interactions be¬ 
tween the DM and SM fields, either through the Higgs 
boson or through direct couplings to quarks and leptons. 
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Fermionic DM {N = 6) 


Scalar DM {N = 6) 



G(l)v 

G(l)s 

G(1)l 

-^,1' 

CSI 

Gdm 

1 

0 

1 

0 

/ 

/ 

UD'^X 

2 

-1/2 

-1 

0 

X 

X 

Q^X 

2 

1/2 

-1 

0 

X 

X 

QD^X 

3 

-1 

0 

-1 

X 

X 

H^LX 

3 

0 

1 

0 

X 

X 

Q'^DX 

3 

1 

1 

0 

/ 

X 

Q^UX 

4 

-1/2 

-1 

0 

X 

X 

Q^X 

4 

3/2 

0 

-3 

X 

X 

L^X 

5 

-1 

0 

-1 

/ 

X 

H^LX 

5 

0 

0 

-1 

/ 

/ 

H^H^LX 

5 

1 

0 

-1 

/ 

X 

H'^^HLX 

5 

2 

0 

-1 

/ 

X 

H^^LX 


Table I. Classification of fermionic DM which is stable up to 
dimension six decays. Columns 1-4 list DM charges under the 
SM symmetry group. Column 5 has a / for DM that remains 
stable up to dimension six decays even when B and L are 
explicitly broken while preserving proton stability. Column 
6 has a / for DM that is safe from spin-independent DM- 
nucleon scattering through the Z boson. Column 7 lists an 
example operator mediating DM decay. 


However, the DM charges are chosen by construction to 
eliminate renormalizable decays, so the latter are highly 
suppressed. On the other hand, the former are allowed 
for scalar DM, and from an effective theory perspective 
should be included if they are permitted by the symme¬ 
tries of the theory. 

If gauge interactions dominate DM annihilations, then 
the relic abundance is fixed by thermal freeze-out W- 
Since the strength of the SU{2 )gauge interac¬ 
tions are known, there is only one free parameter avail¬ 
able to tune the relic abundance: the DM mass. As 
discussed in Hi, to obtain DM of the observed relic 
abundance today requires M of order the weak scale or 
slightly above, depending on the DM spin and quantum 
numbers. Of course, direct couplings to the Higgs US 
can strongly affect the relic abundance. However, anni¬ 
hilations mediated by the Higgs generally add to the DM 
annihilation cross-section, depleting the DM abundance 
and thus requiring a larger DM mass. 

Another approach to DM generation is asymmetric 
DM, which fits quite naturally into the effective the¬ 
ory picture presented here. In models of asymmetric 
DM there is an approximate DM number X that is con¬ 
served along with B or L. In the early universe, some set 
of interactions, for example higher dimension operators 
or sphaelerons, break some combination of X, B, and 
L, thus sharing particle asymmetries among the quarks, 
leptons, and DM. As we have shown, the majority of ef¬ 
fectively stable DM candidates carry non-zero B or L, 
which implies that the very same operators that mediate 
DM decay can also share the asymmetries in the early 
universe. A similar observation was noted about asym¬ 
metric DM models in [B]. Alternatively, these kinds of 


SU{2)l 

G(l)v 

G(l)s 

G(1)l 


C^SI 

Odm 

1 

0 

0 

-2 

X 

X 

H'^L-^X 

3 

0 

0 

-2 

X 

X 

H^L^X 

5 

0 

0 

-2 

X 

X 

H^L^X 

5 

1 

0 

-2 

X 

X 

H^HL^X 

5 

2 

0 

-2 

X 

X 

H^^L'^X 

6 

1/2 

0 

0 

X 

X 


6 

3/2 

0 

0 

X 

X 

h'^^hx 

6 

5/2 

0 

0 

X 

X 

H^^X 


Table II. Classification of scalar DM which is stable up to 
dimension six decays. Notation from Tab. 

higher dimension operators can also be applied in the 
reverse direction to produce the baryon asymmetry m- 

If the decay operator is in chemical equilibrium in the 
early universe, then the DM and B or L asymmetries 
will be shared according to their charges [S]. As is well- 
known, for efficient sharing this typically requires light 
asymmetric DM of order the ~ GeV scale rather than 
the weak scale. For DM that is a gauge singlet, such 
as the first entry in Tab. [l] this offers a viable model of 
asymmetric DM, provided additional annihilation modes 
to deplete the symmetric abundance. 

On the other hand, models with DM carrying SM 
gauge charges are excluded by LEP pn [IH] for ~ GeV 
scale masses. In this case, asymmetric DM is possible 
only if asymmetry sharing through the decay operator 
is inefficient, so the abundance of DM is less than the 
amount prescribed by chemical equilibrium, thus requir¬ 
ing a larger DM mass. 


III. EXPERIMENTAL CONSTRAINTS 

In this section we survey the bounds on effectively 
stable DM from laboratory and telescope experiments. 
As we will see, bounds from direct detection and proton 
decay have an interesting connection to explicit breaking 
of B and L. 


A. Direction Detection 

Experimental bounds on spin-independent DM- 
nucleon scattering are extremely stringent. In particu¬ 
lar, Dirac or complex scalar DM with a spin-independent 
coupling to the Z boson is excluded by many orders of 
magnitude. Such an interaction arises when Y yf 0, so 
direct detection constraints are trivially evaded by DM 
with vanishing hypercharge. 

In Tabs. and |TT] and Tabs. |HI| and |IV[ the column 
denoted ctsi carries a Z for models which have zero hy¬ 
percharge and are thus safe from direct detection, and 
X otherwise. For the latter, models of hypercharged DM 
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Fermionic DM {N = 7) 


Scalar DM {N = 7) 


SUi2)L 

U{1)y 

U{1 )b 

Uil)L 


0-31 

Obm 

3 

1 

-1 

0 

/ 

X 

~WQfFY 

3 

1 

0 

-3 

X 

X 

HL^X 

4 

-3/2 

0 

-1 

/* 

X 

H^LX 

4 

3/2 

1 

0 

/* 

X 

H^Q^UX 

4 

-1/2 

1 

0 

/* 

X 

HQ'^DX 

5 

-1 

-1 

0 

/* 

X 

HQ^X 

5 

0 

-1 

0 

/* 

/ 

H^Q^X 

5 

1 

0 

-3 

/* 

X 

HL^X 

5 

2 

0 

-3 

/* 

X 

H^L^X 

6 

-3/2 

0 

-1 

/ 

X 

H*LX 

6 

-1/2 

0 

-1 

/ 

X 

H^H^LX 

6 

1/2 

0 

-1 

/ 

X 

H^^H^LX 

6 

3/2 

0 

-1 

/ 

X 

H^^HLX 

6 

5/2 

0 

-1 

/ 

X 

H^^LX 


Table III. Classification of fermionic DM which is stable up 
to dimension seven decays. Notation from Tab. [Ij but with 
a /* to indicate models which are cosmologically stable at 
dimension five or less but may still decay at dimension six 
when B and L are violated consistent with proton stability. 


SU{2)l 

t/(l)v 

G(1)b 

U{1)l 

p,]jj 

USI 

Obm 

1 

0 

-1 

-1 

/* 

/ 

Q^LX 

2 

-1/2 

1 

-1 

X 

X 

D^LX 

2 

-1/2 

0 

-2 

X 

X 

H^L^X 

2 

1/2 

1 

-1 

X 

X 

UD^LX 

3 

-1 

-1 

-1 

X 

X 

QU'^LX 

3 

-1 

1 

1 

X 

X 

QD^LX 

3 

0 

-1 

-1 

/* 

/ 

Q^LX 

4 

-3/2 

-1 

1 

X 

X 

Q^EX 

4 

-1/2 

-1 

1 

X 

X 

Q'^ULX 

4 

-1/2 

0 

-2 

X 

X 

H^L^X 

5 

0 

-1 

-1 

/* 

/ 

Q^LX 

5 

2 

0 

-4 

X 

X 

L^X 

6 

-1/2 

0 

-2 

/* 

X 

H^L^X 

6 

1/2 

0 

-2 

/* 

X 

H'^H'^L^X 

6 

3/2 

0 

-2 

/* 

X 

H^^HL^X 

6 

5/2 

0 

-2 

/* 

X 


7 

0 

0 

0 

/ 

/ 


7 

1 

0 

0 

/ 

X 


7 

2 

0 

0 

/ 

X 


7 

3 

0 

0 

/ 

X 

H^^X 


can evade these direct detection if we allow for additional 
structures which we now discuss. 

In particular, limits on spin-independent scattering via 
the Z boson are null if there is an even tiny mass split¬ 
ting S between the components of the Dirac fermion or 
complex scalar [HIID]. In this case, scattering through 
the Z boson is inelastic, requiring additional energy to 
excite the incoming DM particle into a neighboring mass 
eigenvalue. In particular, scattering is kinematically for¬ 
bidden provided the mass splitting exceeds 


where (3 ~ 200 km/sec is the DM velocity and /i is the 
reduced mass of the DM-nucleus system. For the typ¬ 
ical atomic weights of targets such in experiments like 
CDMS HD, XENONIOO |12], and LUX [23], for M - 100 
GeV this translates into a bound 5 ^ 10 keV. For very 
small M, the reduced mass will decrease and so too will 
the lower bound on S, but these regions of light DM are 
excluded by LEP for F 0 [Ul HH] . 

However, inducing the requisite mass splitting 6 re¬ 
quires new operators that explicitly break the DM par¬ 
ticle number associated with Dirac fermions or complex 
scalars. For DM candidates that carry B or L, this im¬ 
plies explicit breaking of B or L. However, the spurion 
responsible for explicit breaking, (AH, AL), enters with 
twice the B and L charge of the DM. Consequently, even 
with these splitting operators, there is still an unbroken 
Z 2 subgroup of B or L that maintains DM stability. 

For example, consider a fermionic DM particle that is 
a doublet of hypercharge Y = 1/2. The leading opera¬ 
tor that can split its components is A mass 


Table IV. Classification of scalar DM which is stable up to 
dimension seven decays. Notation from Tab. [^ 

splitting sufficient to evade direct detection requires 

A < — - 10® GeV, (9) 

0 

which is a low cutoff in the context of DM stability. For 
even larger hypercharges, the level splitting operator in¬ 
volves more Higgs fields, and the requisite cutoff is even 
lower, dropping to ~ 30 TeV for Y = 1. 

According to the philosophy of effective theory, the 
bound in Eq. (|^ defines a cutoff scale at which all 
higher dimension operators allowed by symmetry should 
be present—including those which can mediate DM de¬ 
cay. For dimension five, six, and seven operators, such 
a low cutoff is inconsistent with cosmological limits on 
weak scale DM. This is depicted in Fig. where the blue 
shaded region shows that below above A > 10® GeV, the 
cutoff is too high to induce a sufficient mass splitting 
to evade direct detection, so F = 0 is forbidden. Gon- 
versely, effectively stable DM with non-zero hypercharge 
can only evade direct detection if there is a low cutoff, in 
which case cosmological stability requires the DM decay 
operator be dimension eight or higher. 

In general, if the mass splitting requires a higher di¬ 
mension operator then direct detection is inconsistent 
with the criterion of cosmologically stable DM. On the 
other hand, there is no issue if the splitting operator is 
renormalizable. For example, this is possible for complex 
scalar DM with F = 1/2, which permits the renormal¬ 
izable operator For M = 100 GeV, a splitting 

of (5 = 10 keV is achieved for A > 10“®, which is eas¬ 
ily satisfied for order one couplings. If A carries B or 
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L, this mass splitting operator explicitly breaks B or L 
down to a discrete baryon or lepton parity, however still 
maintaining DM stability. Alternatively, this interaction 
is symmetry preserving if X does not carry B or L, as 
is the case for fc-plets of SU{2)l with large values of k 
which couple only to Higgs bosons in the decay operator. 

Even for models which evade bounds on Z-mediated 
scattering, direct detection may still impose constraints 
on Higgs boson exchange. For example, scalar DM- 
nucleon scattering can occur via the Higgs portal inter¬ 
action, WHX^X |211. For fermionic DM, the analogous 
interactions are higher dimension. Non-singlet fermionic 
and scalar candidates can also scatter with nucleons at 
loop level via multiple gauge boson exchange, which may 
be observable in the next generation of direct detection 
experiments |25j . 


B. Indirect Detection 

Since these DM candidates eventually decay on cosmo¬ 
logical time scales, they are naturally probed by cosmic 
ray telescopes. Conveniently, the authors of [5] studied 
indirect detection constraints on DM decay via high di¬ 
mension operators of the very type considered in this 
paper, albeit with underlying motivation of asymmetric 
DM. In particular, they considered bounds from FERMI, 
PAMELA, AMS-02, and HESS on high energy gamma 
rays and charge particle cosmic rays from electrons, pro¬ 
tons, and anti-protons, obtaining a limit on the DM life¬ 
time which we have taken as a loose input for Eq. 

We refer the reader to for precise numerical bounds, 
but we summarize the salient takeaways below. 

In general, DM carrying B will decay to quarks, yield¬ 
ing anti-protons, while DM carrying L will decay to lep¬ 
tons, yielding positrons and neutrinos. All of the DM 
decays considered here will produce high energy gamma 
rays from charged particle bremsstrahlung, hadronic de¬ 
cays, and inverse Compton scattering of CMB photons 
and starlight. However, due to the large B and L charges 
of most our DM candidates, gamma ray lines are not typ¬ 
ically expected among the theories considered here. An 
exception occurs for DM particles which decay through 
operators involving only the Higgs boson, in which case 
mixing together with a loop of SM particles will induce 
two-body decays of DM to photons. Another possibility 
occurs for DM with unit L number, which can decay to 
photon plus neutrino. 


already places significant constraints on the simplest 
GUTs. From an effective theory viewpoint, proton decay 
is mediated by dimension six operators of the form 


Q^L 
A2 ’ 


( 10 ) 


which breaks B and L but preserves B — L. Current 
limits on proton decay imply a lower bound on the cutoff 
of approximately A > 10^^ GeV [30] i so B and L are very 
well-preserved symmetries. 

As noted earlier, we can parameterize B and L vio¬ 
lation with a dimensionless spurion characterizing hard 
breaking by some number of units of XB and XL. How¬ 
ever, {XB,XL) = (1,1) is of particular note because it 
permits proton decay via the operator in Eq. (10). As 
a result, if we wish to forbid this then we should avoid 
breaking by units of {XB,XL) = (1,1). Remarkably, 
even with {XB,XL) yf (1,1) there are still viable mod¬ 
els of effectively stable DM. In these models there are low 
dimension DM decay operators allowed by the SM gauge 
symmetry, but these operators require explicit breaking 
by {XB, XL) = (1,1) which would also decay the proton. 

For example, this happens in the model described by 
the first row of Tab. [^ Since A is a gauge singlet, the 
gauge invariant operator HLX, if present, would induce 
catastrophically prompt DM decay. However, the exis¬ 
tence of both UD^X together with HLX would require 
symmetry breaking by (XB,XL) — (1,1) which would 
in turn induce dimension six proton decay. Conversely, 
if B and L are explicitly broken while forbidding proton 
decay, an accidental symmetry remains which effectively 
stabilizes DM. 

Note that proton decay is also mediated by oper¬ 
ators beyond dimension six requiring B and L viola¬ 
tion by one unit and an odd number of units, respec¬ 
tively. For example, D^L induces proton decay with 
{XB,XL) — (—1,1) breaking. Forbidding this larger 
class of spurions would allow for even more viable candi¬ 
dates for cosmologically stable DM, but we do not con¬ 
sider this possibility here. 

In Tabs. |l| and |H| and Tabs. |HI| and |IV[ the col¬ 
umn labelled carries a / if DM is still effectively 
stable—that is, cannot decay at dimension five or less— 
even after including any hard breaking spurion with 
{XB,XL) yf (Ij 1)) and X otherwise. The /* entries 
in Tabs. |III| and |IV| indicate models with DM that, while 
stable at dimension five or less, still decays at dimension 
six. For these models, DM is still cosmologically stable 
if the cutoff is higher than was required for dimension 
seven decays. 


C. Proton Decay 

The non-observation of proton decay offers a strong 
motivation for at least approximate B and L conser¬ 
vation. Current limits on —> e“'"7r° and related de¬ 

cay modes from the Super-Kamiokande experiment re¬ 
quire a lifetime of at least ^ 10^^ years |351l29j . which 


IV. SUMMARY AND OUTLOOK 

DM phenomenology often hinges on the assumption 
of a stabilizing symmetry. Naturally, this leads one to 
wonder about the underlying reason for cosmological sta¬ 
bility. In this paper we present an alternative hypothesis 
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whereby DM is long-lived as an accident of the SM sym¬ 
metry group. We have classified all models in which DM 
decay is cosmologically slow and induced by dimension 
six or seven operators. In such cases a sub-Planckian 
cutoff is sufficient to prevent cosmologically prompt DM 
decays which are excluded by indirect detection. All 
candidates for effectively stable DM either carry B or 
L, or reside in a high representation of SU{2)l- We 
have identified those models which are consistent with 
stringent bounds on spin-independent DM-nucleon scat¬ 
tering. Finally, we have accounted for the possibility 
explicit breaking of B and L by arbitrary units. As long 
as the symmetry breaking spurion still forbids proton 
decay, models of effectively stable DM persist. 

Our analysis leaves a number of avenues for future 
work. First, there is the question of DM stabilization 


by SM quark and lepton flavor symmetries, which are 
well-preserved for the light generations. Second, there 
is the question of building these models in a supersym¬ 
metric context. Lastly, it would be interesting to see if 
any of the models presented here arise explicitly in GUT 
constructions. 


ACKNOWLEDGMENTS 

CC is supported by a DOE Early Career Award DE- 
SC0010255 and a Sloan Research Fellowship. DS is sup¬ 
ported in part by U.S. Department of Energy grant DE- 
FG02-92ER40701 and by the Gordon and Betty Moore 
Foundation through Grant No. 776 to the Caltech Moore 
Center for Theoretical Cosmology and Physics. 


[1] A. Arvanitaki, S. Dimopoulos, S. Dubovsky, P. W. Gra¬ 

[15] C. Cheung and D. Sanford, JCAP 1402, 011 (2014) 

ham, R. Harnik, et al., Phys.Rev. D79, 105022 (2009) 

arXiv: 1311.5896 [hep-ph] 

arXiv:0812.2075 [hep-ph] 

[16] C. Cheung and K. Ishiwata, Phys.Rev. D88, 017901 

[2] M. Cirelli, N. Fornengo, and A. Strumia, Nucl.Phys. 

(2013), arXiv: 1304.0468 [hep-ph] 

B753, 178 (2006), arXiv:hep-ph/0512090 [hep-ph] 

[3] M. Cirelli, A. Strumia, and M. Tamburini, Nucl.Phys. 

[17] A. Heister et al. (ALEPH), Phys.Lett. B533, 223 (2002) 
arXiv:hep-ex/0203020 [hep-ex] 

B787, 152 (2007), arXiv:0706.4071 [hep-ph] 

[18] J. Abdallah et al. (DELPHI), Eur.Phys.J. C31, 421 

[4] L. Di Luzio, R. Grober, J. F. Kamenik, and M. Nardec- 

(2003), arXiv:hep-ex/0311019 [hep-ex] 

chia, (2015), arXiv: 1504.00359 [hep-ph] 

[19] L. J. Hall, T. Moroi, and H. Murayama, Phys.Lett. 

[5] D. E. Kaplan, M. A. Luty, and K. M. Zurek, Phys.Rev. 

B424, 305 (1998), arXiv:hep-ph/9712515 [hep-ph] 

D79, 115016 (2009), arXiv:0901.4117 [hep-ph] 

]20] D. Tucker-Smith and N. Weiner, Phys.Rev. D64, 043502 

[6] Y. Zhao and K. M. Zurek, JHEP 1407, 017 (2014) 

(2001), arXiv:hep-ph/0101138 [hep-ph] 

arXiv: 1401.7664 [hep-ph] 

]21] R. Agnese et al. (CDMS), Phys.Rev.Lett. Ill, 251301 

[7] L. Lavoura, S. Morisi, and J. Valle, JHEP 1302, 118 

(2013), arXiv: 1304.4279 [hep-ex] 

(2013), arXiv:1205.3442 [hep-ph] 

]22] E. Aprile et al. (XENONIOO), Phys.Rev.Lett. 109, 

[8] O. Antipin, M. Redi, A. Strumia, and E. Vigiani, 
(2015), arXiv: 1503.08749 [hep-ph] 

181301 (2012), arXiv: 1207.5988 [astro-ph.CO] 

[23] D. Akerib et al. (LUX), Phys.Rev.Lett. 112, 091303 

[9] T. Banks and N. Seiberg, Phys.Rev. D83, 084019 (2011) 
arXiv:1011.5120 [hep-th] 

(2014), arXiv: 1310.8214 [astro-ph.CO] 

[24] R. Schabinger and J. D. Wells, Phys.Rev. D72, 093007 

[10] A. Ibarra, A. S. Lamperstorfer, and J. Silk, Phys.Rev. 
D89, 063539 (2014), arXiv: 1309.2570 [hep-ph] 

(2005), arXiv:hep-ph/0509209 [hep-ph] 

[25] J. M. Cline, K. Kainulainen, P. Scott, and C. Weniger, 

]11] S. Ando and K. Ishiwata, (2015), arXiv: 1502.02007 

Phys.Rev. D88, 055025 (2013), arXiv: 1306.4710 [hep- 

[astro-ph.CO] 

[12] G. Giesen, M. Boudaud, Y. Genolini, V. Poulin, 

ph] 

[26] K. Abe et al. (Super-Kamiokande), Phys.Rev. D90, 

M. Cirelli, et al., (2015), arXiv: 1504.04276 [astro- 

072005 (2014), arXiv: 1408.1195 [hep-ex] 

ph.HE] 

]13] C. Rott, K. Kohri, and S. C. Park, (2014), 

[27] V. Takhistov et al. (Super-Kamiokande), Phys.Rev.Lett. 
113, 101801 (2014), arXiv:1409.1947 [hep-ex] 

arXiv: 1408.4575 [hep-ph] 

[28] K. Abe et al. (Super-Kamiokande), Phys.Rev.Lett. 113, 

[14] J. M. Cline and P. Scott, JCAP 1303, 044 (2013) 

121802 (2014), arXiv: 1305.4391 [hep-ex] 

arXiv:1301.5908 [astro-ph.CO] 

]29] C. Regis et al. (Super-Kamiokande), Phys.Rev. D86, 

012006 (2012), arXiv: 1205.6538 [hep-ex] 

]30] P. Nath and P. Fileviez Perez, Phys.Rept. 441, 191 
(2007), arXiv:hep-ph/0601023 [hep-ph] 




